Using Three-Body Recombination to Extract Electron Temperatures of Ultracold 

Plasmas 
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Three-body recombination, an important collisional process in plasmas, increases dramatically 
at low electron temperatures, with an accepted scaling of T^T 9 '' 2 . We measure three-body recom- 
bination in an ultracold neutral xenon plasma by detecting recombination-created Rydberg atoms 
using a microwave-ionization technique. With the accepted theory (expected to be applicable for 
weakly-coupled plasmas) and our measured rates we extract the plasma temperatures, which are 
in reasonable agreement with previous measurements early in the plasma lifetime. The resulting 
electron temperatures indicate that the plasma continues to cool to temperatures below 1 K. 

PACS numbers: 



Three-body recombination (e - + e~ + A + — > e~ + A*) 
is a fundamental collisional process in plasmas that is 
dominant at sufficiently low electron temperatures due 
to its T(T 9 ^ 2 dependence. In ultracold plasmas (UCPs), 
the observation of copious Rydberg atom production [l[ 
and the observation of T c almost independent of initial 
energies Q show that three-body recombination (3BR) 
and its associated heating play a critical role in the evo- 
lution of UCPs. Early-time T c measurements [2[ and 
simulations 3] suggest that the electrons remain weakly 
coupled, so that traditional 3BR theory is still valid (in 
the strong-coupling regime, the 3BR rate is predicted to 
be reduced below the 9/2 scaling law to a T^ 1 rule 0]). 
A measurement of 3BR in an UCP can thus be used 
to test 3BR theory by using existing T c measurements. 
This is less than ideal, given the paucity of measurements 
and the sensitivity of the rate constant to T D due to the 
9/2 power. Conversely, using 3BR theory, T c can be ex- 
tracted from the measured 3BR rate. This is relatively 
insensitive to the value of the rate constant (due to a 2/9 
power law), and can be used to make the first measure- 
ment of T c throughout the life of the plasma. In addition, 
modifications of the rate constant due to strong-coupling 
will overestimate T , i.e. our extracted T c are an upper 
limit. We note that in addition to furthering our un- 
derstanding of UCPs, a study of 3BR may aid in using 
plasmas with similar parameters (albeit at high magnetic 
fields) to optimize production of anti- hydrogen 

In this work, we measure the instantaneous Rydberg 
atom production rate in an expanding ultracold xenon 
plasma as a function of the time elapsed after plasma 
formation. By applying a short microwave pulse to 
the UCP, the Rydberg population for principal quantum 
numbers N>35 is ionized and subsequently detected by 
a micro-channel plate detector. Using two such pulses 
separated by up to a few microseconds, we measure the 
refilling of the depleted Rydberg levels that form dur- 
ing the interval between the two pulses. In this man- 
ner we determine the instantaneous Rydberg formation 



rate during the plasma expansion. There are several 
processes that may contribute to this refill rate, in- 
cluding 3BR, radiative recombination, blackbody-driven 
transitions from low Rydberg levels to higher (and thus 
microwave-ionizable) levels, and electron-Rydberg colli- 
sions that drive transitions between N- levels. Of these 
processes, 3BR has the strongest inverse-temperature de- 
pendence and will be the most prevalent process at low 
temperatures. Simulations based on the accepted expres- 
sions for these various processes and our T c support our 
assumption that 3BR is the only significant process, so 
that the refill rates we observe are due to 3BR. We mea- 
sure the total 3BR rate, which should be proportional to 

9/2 

K^brT c , where Kssr, is the three-body rate constant. 
We use the accepted K^br to extract T e at later times 
in the plasma, yielding the first T c measurements in an 
UCP later than 10 fis after the plasma creation. 

Our creation of the UCP is similar to previous work 
Q. A magneto-optic trap is used to collect ~ 5 x 10 6 
metastable xenon atoms and to cool them to a temper- 
ature of approximately 20 /iK. The density distribution 
is roughly Gaussian with a rms radius cfq *~ 260 /im and 
a peak density of about uq — 1 x 10 16 m~ 3 . The plasma 
is then produced using a two-photon excitation process 
(882 nm + 514 nm, 10ns pulse), ionizing ~ 30% of the 
atoms. The initial electron energy AE is controlled by 
the frequency of the 514-nm photon; for this work we 
use AE = 3 K. Using a micro-channel plate detector and 
a series of biased voltage grids (the plasma is in a ~ 5 
mV/cm electric field), we detect the emission of electrons 
from the plasma. Figure 1 is an example of such a signal. 

The plasma is unconfined, expands freely into vacuum, 
and lives for ~ 200 fis. As it expands, it maintains a 
roughly Gaussian density distribution (confirmed by ion 
imaging, to be published); we assume that 
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n ~2x 10 15 m" 3 @, H, HI- While simulations indicate 
some deviation from this distribution fiol . we note 
that this deviation should have only a small effect on our 
application of cq. 1. Furthermore, although the plasma 
is constantly losing both ions and electrons, n D and «i 
will remain approximately equal, especially at the plasma 
center. Most deviations from the Gaussian distribution 
occur in the outer, low-density region, but since almost 
all 3BR occurs in the high-density central region of the 
plasma due to the n^rii dependence of the 3BR rate: 

= K 3BR T- Q / 2 J ' nl(r) ni {r)Aiir 2 dr, (2) 

eq. 1 is a good approximation for use in 3BR calculations. 

In addition to this t- and r-dependent density distri- 
bution, the plasma also exhibits time-dependent elec- 
tron and ion temperatures due to heating effects (3BR, 
disorder-induced heating) and cooling effects (adiabatic 
expansion, electron evaporation). Effective ion temper- 
atures (distinct from thermal temperatures) have been 
directly measured 0, EH , but only during the first ~ 4 
/is of plasma expansion. While the thermal electron tem- 
perature, T e , has been estimated using simulation results 
[id HH , measurements have only been made for the first 
10 (is of plasma expansion Q. The use of Tonks-Dattner 
resonances as a plasma diagnostic may provide a way to 
measure T c at times up to ~ 35 fxs 13]; however, cur- 
rent theory prevents it from being a quantitative mea- 
surement. Our 3BR measurements span a much larger 
region of the plasma lifetime than T e measurements, al- 
lowing only limited comparison of our measured rates to 
calculated 3BR rates using known T c . We note that sig- 
nificant Rydberg production occurs throughout the entire 
lifetime of the plasma [l| , allowing temperature informa- 
tion to be determined throughout. 

Microwave ionization of Rydbergs is a non-resonant 
process and largely independent of the microwave fre- 
quency; we use 2.4 GHz pulses to take advantage of an 
electrical resonance in our apparatus. We find that these 
pulses ionize Rydbergs with N>35 by applying a DC field 
ionization ramp immediately after a microwave pulse and 
noting that the signal corresponding to N>35 Rydbergs 
is almost entirely eliminated. We estimate a microwave 
field amplitude of ~ 220 V/cm by observing which prin- 
ciple quantum numbers are ionized, as we are unable to 
directly measure the microwave field at the plasma. We 
find that we are not directly heating the plasma elec- 
trons. Starting with a plasma electron energy of AE= 
300 K (to suppress Rydberg formation) and applying the 
microwave pulses at varying times in the plasma, we ob- 
serve that the plasma expansion and the electron emis- 
sion signal do not substantially change with the addition 
of the microwave pulses. The microwave pulses can heat 
the plasma indirectly through the ionization electrons, 
as those electrons undergo a few collisions while exiting 
the plasma. The electrons resulting from microwave ion- 
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FIG. 1: Electron emission signal from an expanding UCP, 
averaged over 40 runs. The double peaks at ~ 55/xs are the 
response due to a pair of short (100 ns) microwave pulses. 
The first of the two pulses is held fixed and ionizes Rydbergs 
that have formed in the plasma. Following this pulse, the 
Rydberg population is refilled by three-body recombination. 
The second pulse is applied at varying times and the num- 
ber of Rydbergs ionized is counted as a function of the time 
between the two pulses (inset). 



ization are warm (>10 K) [3[ with respect to the plasma 
electrons, but leave the plasma quickly, with limited ther- 
malization with the plasma electrons. We minimize this 
effect by keeping the microwave pulse duration as short 
as possible (100 ns) while still fully ionizing the Rydberg 
population. 

We apply the first microwave pulse at a chosen time 
in the plasma expansion in order to remove the existing 
N>35 population. We then apply a second pulse after a 
variable delay. This results in a double-peaked electron 
signal (fig. 1 near 55 /xs). We fit to the broad emis- 
sion curve and use this as a background to determine the 
number of electrons contained in the second peak[l4T|. A 
single microwave pulse (as well as DC field ionization) 
can also provide a means to measure the number of Ryd- 
berg atoms in the plasma, but the long Rydberg lifetimes 
and the collisional redistribution of Rydberg states make 
this measurement difficult to interpret in the context of 
3BR, as it is the result of past Rydberg production and 
redistribution. 

We note that the width of a electron emission peak due 
to our short microwave pulse is a few times wider than 
the duration of the microwave pulse. There is also a delay 
before we begin to detect Rydberg refilling (fig. 1 inset). 
If the plasma is heated a small amount, as is likely follow- 



ing a microwave ionizing pulse, 3BR would be inhibited 
immediately after an ionizing pulse, explaining the de- 
lay before we observe Rydberg formation. This is not a 
major effect, however, because the energy put into the 
plasma electrons by the newly-ionized electrons should 
largely be removed through the loss of hot electrons (as 
observed in our electron emission peaks). The increased 
width of the electron emission peaks is likely the result 
of a Coulomb-driven spread of the electron bunch during 
the time of flight to the detector. 

Using this two-pulse method, we obtain the Rydberg 
atom refill rate throughout the plasma expansion. A typ- 
ical refill curve (for t = 55 /is) is shown in the inset plot 
of figure 1. We fit a straight line to the first few hun- 
dred nanoseconds of each refill curve and use the slope 
of that line as the short-time refill rate (that is, the in- 
stantaneous 3BR rate). At late times (> 1-2 zxs) in the 
refill curve the rate decreases, likely due to the plasma 
regaining a near-equilibrium Rydberg distribution. 

Figure 2 shows the measured Rydberg refill rates. 
Because multiple processes can create Rydbergs in the 
N>35 states, we considered the possibility that these 
measured rates may be dominated by non-3BR events. 
Using established rates for radiative recombination [l5j |. 
blackbody-driven transitions and ionization [l6j |. radia- 
tive decay 0, and 3BR 0, we simulate [Hj] the change 
in the populations in each TV level for 3<N<200 for sev- 
eral [is starting at each time we obtained refill rates, using 
the plasma densities at that time (eq. 1) with varied T c . 
We use an initial Rydberg population distribution which 
is truncated at N=35 to simulate the microwave ioniza- 
tion event, then include all of the above-cited processes 
in a Monte-Carlo model of the time-dependent Rydberg 
distribution. The results of this simulation confirm that 
3BR is the process most responsible for creating Ryd- 
bergs in N>35 states following an ionization pulse, as 
the other processes together contribute less than a few 
percent to the rates we measure. The simulation also 
supports our measurement of only N>35 Rydbergs as a 
complete 3BR measurement, as 3BR at such low tem- 
peratures almost exclusively populates high levels (N> 
80). 

If we then assume that 3BR theory with K-^-qk = 
4.5 x 10~ 21 m 6 K 9 / 2 s -1 provides us with an accurate rate 
equation (eq. 2) [13], we can use our measured recom- 
bination rates together with the Gaussian expansion ap- 
proximation (eq. 1) to extract T e (t), plotted in fig. 3a. 
We note that these T c (t) are similar to those predicted by 
simulations [H, [13] • A simple power-law fit of the temper- 
ature curve of figure 3a gives T c (t) ~ i _12 W. Adiabatic 
cooling would suggest a T e (t) ~ t~ 2 relationship; the dif- 
ference is likely due to the significant heating effects from 
3BR. This method of determining the electron tempera- 
ture finds that the electrons are well below one Kelvin at 
later times, by far the lowest T c observed in an ultracold 
plasma. Indeed, extrapolating to the end of the plasma 
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FIG. 2: The total Rydberg refill rates. The early increase 
in the rate is likely due to 3BR turning on as temperatures 
drop; the rate peaks and then drops to a constant value as the 
density and temperature continue to decrease. Note that the 
peak plasma density falls by over two orders of magnitude, 
from lxlO 14 m~ 3 at 10 fis to 6x10" m~ 3 at 60 iis, while 
the 3BR rate only decreases by a factor of two. The error 
bars represent the la standard uncertainty resulting from the 
linear fits to the refill curves (such as the inset of fig. 1). 

lifetime (~ 200 /xs), the power-law fit would indicate tem- 
peratures as low as ~ 200 mK. While we observe Rydberg 
production throughout the lifetime of the plasma, we did 
not measure past t=65 /xs because the indirect heating 
effects of the 100 ns microwave pulses become significant 
compared to T e , interfering with extraction of the 3BR 
rates. We note that T e is rather insensitive to the value 
of -K3BR, as a factor of two change in the rate constant 
results in only a 16% change in temperature. 

The temperature in figure 3a is for a plasma with an 
initial T c of 3 K. Also plotted are previous T e measure- 
ments for an initial AE=10 K (squares, |2|) and simula- 
tion results for a AE=66 K (triangles, [3j). Despite the 
different initial temperatures, the three sets are roughly 
consistent with one another; as seen in the T e tend 
to converge to similar values despite very different AE's. 

Using the temperature calculations of figure 3a with 
the density expression given in eq. 1, we calculate the 
electron coupling parameter A = e 2 /(47reoa/cfcT e ), where 
a is the average Wigner-Seitz length. A(i) is plotted in 
figure 3b. The system remains weakly coupled (A -C 1), 
but the coupling is increasing as the plasma expands and 
cools. 

Some work has been done to modify eq. 2. Using 
a Michie-King distribution for the electron velocity dis- 
tributions instead of a Maxwell-Boltzmann distribution 
results in an average corrective factor to eq. 2 of 1 to 1.1 
[HI, increasing our temperatures less than 2.5 %. Mod- 
ifications to eq. 2 to correct for low temperatures by 
using a density dependent cutoff for the Rydberg levels 
created by 3BR |4( do not apply to our plasma. The the- 
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ory modifications are only applicable for a plasma denser 
than ours by several orders of magnitude (our plasma is 
still in the weakly-coupled regime) . Any strong-coupling 
correction is likely to result in lower calculated 3BR rates 
at a particular T c . This indicates that if such a correc- 
tion to 3BR rates is needed for the parameters of this 
UCP, our calculated T c would be an upper bound esti- 
mate of the actual plasma T c . It should be noted that 
even a substantial change to the 3BR rate equation could 
be masked in the temperature plot of figure 3a due to the 
exponent of 9/2 on the T c term. 

Recent work 20] indicates that at electron tempera- 
tures below ~ 1000 K, quantum 3BR should populate 
lower N levels than classical 3BR theory would predict. 
Classical 3BR theory predicts almost no formation in 
N<35 states. If quantum 3BR effects result in the for- 
mation of N<35 Rydbergs, the rates we measure in this 
work would undercount the actual 3BR rate. We can 
not directly compare our results with the calculations 
presented in [2(J, as those calculations did not go lower 
than T c ~ 100 K, nor have we made comprehensive mea- 
surements of the N-level distributions that are created in 
the 3BR-driven refilling process. However, our apparent 
agreement with classical 3BR calculations suggests that 
quantum 3BR at T ~ 1-10 K has little effect on the 
cumulative N>35 population. 

We have directly observed the Rydberg formation rate 
at varying times in the expansion of an ultracold plasma 
by using a short microwave pulse to ionize a Rydberg 
population, followed by a second pulse to probe the re- 
filling of that population. Simulations show that 3BR is 
the mechanism by which this Rydberg population refills, 
because other mechanisms are suppressed due to the low 
temperatures in an ultracold plasma. In future work, we 
can measure the N-level distribution that is formed by 
3BR to determine the importance of quantum effects on 
3BR theory [13] , by varying the microwave field strength 
of the second pulse. 

Although standard 3BR theory is expected to fail in 
strongly coupled systems, we use a typical 3BR expres- 
sion and extract T c (t) for the plasma. Because the result- 
ing T e estimates match well with another measurement 
of the temperature and with simulation results (at least 
at early times in the plasma expansion), the use of the 
accepted 3BR rate expression in this work appears justi- 
fied. Our method is self-consistent in that the extracted 
temperatures show the plasma to still be weakly coupled. 
Future measures of T c in an independent manner would 
allow a direct test the validity of the 3BR expression in 
UCPs. This is the lowest-T e measurement of three-body 
recombination rates to date; the measured rates indicate 
the plasma system achieves sub-Kelvin electron temper- 
atures at late times. 

This work was partially supported by the National Sci- 
ence Foundation PHY0245023. 
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FIG. 3: (a) We calculate T c using established 3BR theory 
with the rates we measure and an assumed self-similar Gaus- 
sian expansion of the plasma (diamonds; also inset). The 
initial energy is AE=3 K. For comparison, we plot earlier 
T e measurements for AE=10 K (squares) [2j and simulation 
results for AE=66 K (triangles) (b) The corresponding 
plasma coupling parameter A. The error bars are the propa- 
gated uncertainties from figure 2. 
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